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ON  RECENT  PROGRESS  IN  AND  UTILIZATION  OF  ASTROGEODETIC- INERTIAL  AND  ASTROGEO- 
DET I C -GRAD I OMETR I C  GRAVITY  VECTOR  DETERMINATION 

ZU  JUNGSTEM  FORTSCHRITT  IN  DER  UND  ANWENDUNG  VON  ASTRQGEODATISCH- INERTIALER 
UNO  ASTROGEQDATI SCH-GRAD IOMETR I SCHER  SCHWEREVEKTORBESTIMMUNG 

PROGRES  RECENTS  ET  UTILISATION  DU  CALCUL  DU  VECTEUR  PESANTEUR  PAR  ASTRO- 
GEODESIE  INERTIELLE  ET  ASTROGEODESIE  GRADIQMETRIQUE 

H.  Baussus  von  Lueczow 


SUMMARY 

Following  an  introductory  status  review,  the  paper  presents  a  new  astrogeo- 
detic-inertial  method  for  vertical  deflection  determination  under  considera¬ 
tion  of  horizontal  channel  interactions  and  suitable  for  the  utilization  of 
single  or  multiple  track  data  in  conjunction  with  supporting  astrogeodetic 
deflections.  It  then  discusses  the  problem  of  optimal  deflection  area  adjust¬ 
ments  under  consideration  of  ocher  work.  Subsequently,  accuracy  aspects  of 
astrogeodetic-gradiometric  gravity  vector  determination  based  on  the  Bell 
gravity  gradiometer,  the  simultaneous  utilization  of  gradiometric  and  inertial 
data  in  terrestrial  surveys,  and  potential  applications  are  addressed. 

ZUSAMMENFASSUNG 

Im  Anschluss  an  eine  einfCfhrende  Statustlbersicht  beschreibt  dieser  Artikel 
eine  neue  astrogeodatisch-inertiale  Mechode  fUr  die  Bescimmung  von  Lotab- 
weichungen  uncer  Beru'cksichtigung  gegensei tiger  horizontaler  Kanaleinfldsse 
und  geeignec  fUr  die  Verwendung  von  Einzel-oder  Mehrf achtrassendaten.  Alsdann 
ero'rtert  er  das  Problem  opcimaler  Fiachenausgleichung  unter  Betrachtung 
anderer  Arbeiten.  Anschliessend  wird  auf  Genaugkeitsaspekte  astrogeoddtisch- 
gradiometrischer  Schwerevektorbestimmung ,  die  gleichzeitige  Verwendung 
gradiometrischer  und  inertialer  Dacen  bei  Messungen  auf  der  Erde  und 
Anwendungsmo'glichkeiten  eingegangen. 

RESUME 

Apres  avoir  rappele  l'historique  de  la  question,  cet  article  presence  uiie 
nouvelle  mechode  d ' ascrog£od£sie  inertielle  pour  le  calcul  de  la  deflexion 
verticale  en  tenant  corapte  des  interactions  des  canaux  horizontaux  et 
permettanc  1 ' uti lisation  des  donn£es  d'une  ou  plusieurs  cransversales 
con joinceraenc  avec  les  deflexions  ascrog£od£siques.  II  examine  ensuite  le 
problfeme  des  meilleurs  r£glages  en  macffere  de  deflexion  pour  d'autres 
emplois.  Enfin  sonc  presentees  1' importance  de  la  precision  pour  la  raesure 
ascrogSodesique  et  gradiometrique  du  vecteur  pesanteur  £  l'aide  de  l'appareii 
3ell  de  mesure  du  gradianc  de  pesanteur,  l'emploi  simultane  des  donn£es 
gradiometriques  et  inertielles  pour  les  releves  copographiques  terrestres 
ainsi  que  d'autres  emplois  possibles. 


611.2/2 


1 .  INTRODUCTION 

The  latest  consolidated  presentations  concerning  gravity  vector  determination 
by  means  of  astrogeodetic-lnertlal  and  astrogeodetlc-gradlometrlc  data  can  be 
found  in  the  proceedings  of  the  Second  International  Symposium  on  Inertial 
Technology  for  Surveying  &  Geodesy,  held  in  Banff,  Canada,  June  1-5,  1981.  At 
chat  time,  Adams  and  Hadfield  (1981)  addressed  GEO-SPIN/ IPS-2  improvements  for 
precision  gravity  mesuremenc*,  Bose  and  Huddle  (1981)  described  a  regional 
adjustment  method  for  gravity  vector  determination,  Harris  (1981)  presented 
IPS-2  test  results,  Todd  (1981)  reviewed  modified  IPS-1  test  data,  Heller 
(1981)  elaborated  on  prospects  for  gradiometric  aiding  of  inertial  survey 
systems,  Meczger  and  Jincltano  (1981)  summarized  the  application  of  Bell's 
gravity  gradlometer  and  gravity  meters  to  airborne  and  land  vehicle  gravity 
surveys,  Trageser  (1981)  gave  a  floated  gravity  gradlometer  status  report2, 
Paik  (1981)  reported  on  the  superconducting  gravity  gradlometer,  and  Brown 
(1981)  dealt  with  methods  of  processing  gradlometer  data  for  geophysical 
applications.  With  respect  to  improved  astrogeodetic-inertial  gravity  vectc 
decerminacion,  Honeywell  (1981)  studied  modified  hardware  and  software 
requirements  to  achieve  high  deflection  change  accuracies  of  the  order  of  0. 
arcsec  average  rms  for  a  survey  length  of  about  60  km  with  present  hardware  .. 
the  concext  of  multiple  traverses  (area  adjustment)  and  of  the  order  of  0.1 
arcsec  rms  following  installation  of  improved  accelerometers  and  velocity 
quantizers.  Litton  (1982)  claimed  equivalent  accuracies.  As  to  an  average 
rms  requirement  of  0.1  arcsec,  Litton  proposed  utilization  of  ^creened  G  300 
G2  gyroscopes  with  correlated  random  noise  parameters  of  2.10"  deg  hr~  and 
min  (correlation  length)  and  A  1000  accelerometers  with  corresponding  random 
parameters  of  1  mgal  and  2  min.  Baussus  von  Luetzow  (1982)  developed  a 
coupled  horizontal  channel  optimal  method  for  vertical  deflection  determina¬ 
tion  in  serai-flac  terrain  in  che  context  of  Litton's  local-level  system.  In 
the  field  of  gravity  gradiomecry,  White  (1980)  researched  error  models  and 
related  aspects,  and  Chan  (1982)  presented  progress  in  che  development  and 
testing  of  the  superconducting  gravity  gradiometer.  The  emphasis  In  this 
paper  is  on  che  analysis  of  a  framework  for  optimal  .astrogeodetic-inertial 
determination  of  vertical  deflections  under  consideration  of  the  area  adjust¬ 
ment  problem  including  che  smoothing  method  designed  by  Bose  and  Huddle 
(1981),  and  on  a  discussion  of  significant  aspects  of  gravity  vector  determi¬ 
nation  under  application  of  che  Bell  system,  simultaneous  utilization  of 
gradiometric  and  inertial  data,  and  potential  applications  under  inclusion  of 
che  superconducting  gravity  gradiometer.  Apart  from  technical  information  of 
interest,  this  presentation  has  been  designed  to  arrive  at  some  useful  conclu¬ 
sions  . 
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? 
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IPS-2  relates  to  Honeywell's  inertial  positioning  system  incorporating  a 
space-stable  pLacform,  IPS-l  stands  for  Litton's  corresponding  local-level 
system. 

Draper  Laboratories  floated  gravity  gradiometer  development  and  Hughes 
rotating  gravity  gradiometer  development  were  suspended  in  1980  and  1979, 
respectively. 

The  parameters  are  approximate.  Acer _ _ 
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2.  COUPLED  HORIZONTAL  CHANNEL  OPTIMAL  DETERMINATION  OF  VERTICAL  DEFLEC¬ 
TIONS  IN  SEMI- FLAT  TERRAIN 


Higher  accuracy  requirements  necessitating  installation  of  gyros  and  acceler¬ 
ometers  with  small  error  variances  and  short  correlation  times  and  of  high 
accuracy  velocity  quantizers,  Identified  in  section  1,  require  the  integration 
of  the  system  of  pertinent  differential  equations  for  best  possible  estimation 
of  deflection  components.  In  this  context,  consideration  of  constant  gyro 
biases  require  ultimately  a  Wiener-Kolmogorov  optimization  since  optimal 
Kalman  smooching  is  only  possible  in  Che  case  of  observable  linear  aggregates 
of  random  errors.  The  pertinent  local-level  equations  are: 


X  * 

Vz  -  **N  +  *n  +  aE  -  aEo 

n) 

7  - 

"Vz  ♦  -  g5  +  aN  -  aNo 

$  * 

z 

_i  2 

R  tne.x  +  R  (.T^  +  o^  sec  p)y 

*S  * 

R-1x 

“Ve  ** 

®E  ' 

-R_1y 

”Vz  ***** 

4-y 

(5 

Symbols  used  above  including  total  time  derivatives  of  first  and  second  order 
are : 


'z 

’N 

Sr 

Sc 

SN 

aE*aEo 


V  “No 


i 


meridian  vertical  defleccion 
prime  vertical  defleccion 
vertical  gravity  vector  component 
earth's  mean  radius 
geographic  lacicude 
earth's  inertial  angular  velocity 
system's  east  velocity 
system's  north  velocity 
azimuth  placform  attitude  error 
platform  tilt  error  about  north  axis 
placform  tilt  error  about  east  axis 
east  acceleration  of  survey  vehicle 
north  acceleration  of  survey  vehicle 

correlated  east  accelerometer  error  under  consideration  of  initial 
calibration 

correlated  north  accelerometer  error  uner  consideration  of  initial 
calibration 

-  2Z  *  r1  Vx  too  vertical  spatial  race,  •  D  sin  o 
»  .~*j  -  VjjR”*  north  spatial  rate,  r.»|  *  .1  cos  o 

*  -R"lV  east  spatial  race 
azimuth  axis  angular  drift  race  error 
north  axis  angular  drift  rate  error 
east  axis  angular  drift  race  error 
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For  land  vehicles,  eqs.  (3)  -  (3)  may  be  simplified  by  omission  of  o^,  uc,  and 
pz  -  R  Vx  tn$,  and  by  use  of  constant  accelerations  S»j  and  SE  which  should  be 
approximately  achieved.  Then,  ui2  ■  nsin$  and  »  ftcosq.  The  initial 
conditions  at  t0  »  0  are,  under  consideration  of  plumbline  levelling,  u2(0)  » 
o,  bN(0)  -  n0,  *E(0)  -  ;0,  x(o)  -  y(o)  ■  xQ  -  y(o)  ■  0.  The  accelerations  SN 
and  SE  can  probably  be  neglected.  If  the  system  Is  treated  as  one  with 
constant  coefficients  as  a  good  approximation,  a  closed  solution  as  a  function 
of  time  is  possible.  Because  of  intermittent  Kalman  filter  corrections  and 
the  need  for  numerical  weight  factors,  appears  to  be  advantageous  to 
attempt  a  numerical  soluclon  for  x  and  y  under  utilization  of  terminal 
deflection  and  azimuth  data.  For  economic  reasons  it  is  assumed  that  the 
survey  vehicle  travels  approximately  at  a  constant  speed  when  in  motion  and 
stops  every  3  minutes  for  1  minute.  The  speed  should  not  exceed  10  msec-  in 
order  to  restrict  the  length  of  travel  intervals. 

Solutions  for  lz,  *,  y,  *,  y  are  obtained  in  accordance  with  the 

integration  schemes' 


with  At  »  30  sec  and  possibly  60  sec.  The  solution  structure  at  the  end  of 
the  first  stop  interval,  indicated  by  the  subscript  s  *  I,  is 

■J  V  7  i  j 

r  ,  »  "  a  v*"bw+rca-*-  d3+'"ey  <’■ 

vl  J  v  v  j-  v  v  l  J  v  -  o  v  v  v 

In  eq.  (8),  v  -  g  ( n  -  nQ)  +  a£  -  aEo,  v  -  -g(E  -  *  a,,  -  aNQ. 

Under  consideration  of  Kalman  filter  tilt  corrections,  assumed  for  simplicity 
to  eliminate  the  integrated  first  two  terms  in  eq.  (3)  and  the  first  term  in 
botn  eqs.  (4)  and  (5)  , 
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3 
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T 
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(9) 


w  tie  re  and  subsequent  Ts  represents  an  aggregate  of  tilt-induced  random 

errors.  During  the  stop  interval,  eqs.  (3)  -  (5)  are  integrated  with  x  ■  y  3 
y  *  j  with  a  resultant  effect  on  eqs.  (I)  and  (2).  The  integration  is  then 
continued  and  yields  under  utilization  of  average  biases  the  solution 
structure 


•  “s  A  7  -  "S  5  w  +  C  3*0  5  +  E  Y  R  *  T 

■iS  -  ,  1  I  *  ,  >  >  >  1  J  VS  S 

7»i  .1*1 


where  R^s  designates  residual,  random-type  terms  including  5  -  i.},  3  -  2^ 


At  the  termination  of  the  survey,  when  *n,  ln,  and  Arj  are  available,  it  is 
possible  to  sol /e  for  the  gyro  biases  in  the  form 
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where  the  F's  are  computable,  the  e's  are  an  aggregate  of  v  and  werrors,  and 
the  o's  are  aggregates  of  gyro  random  and  tilt  errors.  Substitution  of  5,  3, 

7  from  eq.  (11)  in  eq.  (10)  and  separation  of  v  and  w  into  def  lection.and 
accelerometer  errors  leads,  under  restriction  to  solutions  for  x  and  y,  to  the 
final  results 
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The  last  4  terms  on  the  right  side  of  both  eqs.  (12)  and  (13)  are  aggregates 
of  random  errors  associated  with  gyros,  accelerometers,  tilt  corrections,  and 
initial  and  terminal  deflection  errors.  Equations  (12)  and  (13)  are 
reformulated  as 


n  n 

n  +  "  k  n  +  1  -  S  -  M  -N 

s  “  j  o  -  ,  j  j  n  n  n 
v» 0  v»0  s  s 


(14 


n  n 

m  n  +  n  \  «  S,  »  M. 

“  „  J  V  u  .  J  ■ J  %  % 

j*0  j* 0  s 


(15) 


where  S,  M,  N  denote  signal,  measurable  message,  and  non— neasurable  noise, 
respectively. 


Under  utilization  of  vertical  deflection  covariance  functions,  the  prime 
dec  lection  can,  for  example,  be  optimally  estimated  in  the  form 


<  lb 


where  i  *  s,  A^e  and  are  matrices  of  regression  coe f f ic ients ,  and  the 

terms  in  brackets  are  message  matrices.  ifith  k  ■  0,1 . .  A{e  and  8ie 

,nn  be  determined  from  the  equations 
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In  eqs.  (17)  and  (18),  where  the  bar  symbol  stands  for  covariance,  the  noise 
covariances  need  only  be  computed  once. 

Simplified  solutions,  particularly  in  the  case  of  approximately  straight 
traverses,  are  possible.  If  averaged  message-type  data  from  repeated  surveys 
are  employed,  the  instrument-generated  noise  covariances  in  eqs.  (17)  and  (18) 
are  to  be  reduced. 


3.  OPTIMAL  REGIONAL  VERTICAL  DEFLECTION  ADJUSTMENT 


The  single  channel  Wiener-Kolmogorov  solution  of  the  astrogeodetic-inertial 
vertical  deflection  determination  problem  was  first  outlined  by  Baussus  von 
Leutzow  (1981)  and  also  considered  by  Litton  (1982).  It  simultaneously  repre¬ 
sents  the  optimal  framework  for  a  regional  adjustment  under  consideration  of 
mulciple  track  data.  In  this  respect,  it  is  advantageous  to  conduct  surveys 
along  approximately  parallel  traverses  and  cross  traverses,  preferably  with  ■ 
new  system's  calibration  at  the  start  of  a  new  traverse  and  without  signifi¬ 
cant  and  rapid  course  changes  in  order  to  reduce  the  effects  of  platform 
heading  sensitivity.  In  practice,  the  computation  of  regression  coefficients 
for  a  particular  5,  n-solution  by  means  of  the  system  of  equations  (17)  and 
( 18)  requires  only  consideration  of  up  to  50  measurements.  Spatial  signal 
covariance  functions  may  only  be  employed  in  moderately  mountainous  terrain. 


Bose  and  Huddle  (1981)  developed  a  different  regional  adjustment  system  free 
of  the  use  of  empirical  signal  covariance  functions.  In  their  approach,  they 
impose  an  orthogonality  restriction  on  the  covariance  of  Fourier  coefficients 
appearing  in  ?QTe  where  Tg  and  T,  are  disturbing  potentials  at  points  P0  and 
Pg.  This  permits  the  representation  of  TQTe  in  terms  of  eigenfunctions  if 


■  2T  .2X  .  -ifr  ri^T 

o“  e  3zz ' o  3zZ/e 


—  4.  — ‘i  '2c.  4-  2H' 

3x  3v'o  3x  3y 'e 


a2  *  const. 


(14 


Within  a  limited  rectangular  area  containing  gridded  message  data  nj,  dgi 
it  is  then  possible,  under  general  assumption  that  T  vanishes  on  the  boundary 
and  chat  the  error  variances  relating  to  the  foregoing  variables  are  constant 
throughout  Che  whole  domain,  to  establish  weight  factor  solutions  in  the  form 


a  -, 

el  1 
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el  l 


'ei  Ugi 


( 20) 


The  above  meChori  suffers  from  Che  serious  restriction  (19)  which  is  less 
accurate  chan  empirical  covariance  functions  derived  under  the  assumption  of 
huraogene l cy  and  isotropy  concerning  TQTe .  As  to  the  structure  thereof,  refer¬ 
ence  is  made  to  Jordan  (1931).  In  conjunction  herewith,  the  assumption  T  »  0 
on  the  boundary  is  incompatible  with  a  finite  domain.  The  requirements  of 
homogeneous  error  variances  and  of  regularly  spaced  gridded  data  are  also  not 
realistic  or  not  easiLy  Co  be  achieved,  respectively. 

4.  ACCURACY  ASPECTS  OF  THE  BELL  GRAVITY  G RADIOMETER  INSTRUMENT  ( GGI ) 

Beil  Aer  ispai  e  TEXTRON,  under  the  technical  leadership  of  E.  Metzger, 
completed  development  of  Che  ball-hearing  rotating  accelerometer  gravity 
,r  lUiomeCc-r  instrument  (GCI)  in  1981.  Its  operational  principle  is  evident 
:  r  «n  the  til. owing.  Be  1 1- turn  i  shed  figure. 
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FIGURE 

Schematic  Illustration  of  Rotating  Bail  GGI  Fixture 


During  1981  3*1.1.  also  concluded  investigations  concerning  gimbal  and  survey 
venicle  self  gradients.  Preliminary  GGI  tests  on  board  a  ship  have  been 
encouraging.  Aircraft  applications  simultaneously  require  continued  position¬ 
ing  daca  as  a  function  of  time.  In  the  context  of  land  vehicle  applications , 
accurate  positioning  is  facilitated  by  Kalman  filter  error  control  under 
utilization  of  velocity  errors  observed  during  periodic  vehicle  stops,  possi- 
oi/  augmented  by  terminal  position  information.  The  total  GGI,  including 
rac<s.  and  power  ‘iipply,  is  considerably  more  complex,  voluminous  and  heavier 
taan  a  pure  inertial  system  and  accordingly  more  expensive  as  to  acquisition 
ma  operation.  intrinsic  GGI  errors  consist  of  instrument  self-generated 
noise  in  trie  platform  environment  ana  environmental  sensitivities  (accelera¬ 
tion,  pressure,  temperature,  magnetic,  etc.)  Self-noise  power  spectral  densi¬ 
ties  for  t:ie  GGI  mounted  at  the  'imbreila  angle  (spin  axis  *  95°  from  horizon¬ 
tal;  nave  been  approximated  in  the  open  literature  by 

*„<:')  -  —  -  3  (21) 

F2 

where  f  denotes  frequency  In  cycles  per  second,  A  is  a  constant  relating 
to  low  frequency  or  red  noise,  B  is  a  constant  indicating  high  freauency 

"  "  _q  -2  -l 

or  white  noise,  E  Is  Eotovos  unit."  10  sec  ‘ ,  and  Hz  is  Hertz  •  1  sec  . 

White  ( 1980)  listed  A  -  2.10’*  E2.Hz  and  16. 10-6  E2.  Hz  *or 

configurations  with  vertical  and  horizontal  spin  axis,  respectively 

and  corresponding  B  •  81  E.Hz”1  and  86  E  .Hz  ”  .  As  to  environmental 

sensitivities,  linear  acceleration  influences  are  the  most  severe, 

particularly  at  multiples  of  the  nominal  rotation  frequency  2  »  n.25  Hz. 

Environmental  errors  represent  approximately  white  noise  and  thus  add  to  B 

in  equation  (21).  The  low  frequency  error  may  he  essentially  Interpreted 

as  random  drift.  In  the  form  (21),  7  does  not  permit  a  Fourier  cosine 

transform  to  determine  the  associated  covariance  function,  “owever,  in  a 

certain  low  frequency  range,  the  first  term  in  eq.  (21)  may  be  replaced  by 

A  (a+f2)“l  where  a  is  an  appropriately  computed  constant.  This  results  in 


.1 
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an  exponential  covariance  function.  The  total  error  variance,  both  self¬ 
generated  and  environmental,  is  in  excess  of  3E  .  Essentially  different 
moving  GGI  information  becomes  essentially  available  every  2  seconds.  Gra¬ 
dient  filters  are  nominally  fourth  order  Butterworth  filters  to  attenuate  1 
3  .1  and  other  noise  manifestations.  Signal  information  is  obtained  as  a  ten 
second  moving  window  average.  For  a  land  vehicle  with  a  speed  of  10  msec-  , 
usable  messages  are  thus  generated  every  20  meters.  Because  of  nonstationarv 
information  in  strongly  mountainous  areas  it  is,  therefore,  necessary  to 
cravel  at  low  speeds,  i.e.,  at  5  msec”  .  For  the  same  reason,  aircraft  over 
strongly  mountainous  regions  should  travel  at  lowest  possible  speed.  Elimina¬ 
tion  of  systematic  drifts  and  biases  has  to  be  accomplished  at  the  start  of  a 
mission.  In  this  respect.  Initial  gravity  tensor  information  is  of  considera¬ 
ble  value.  Otherwise,  observable  gravity  vector  component  errors  at  mission 
termination  would  consist  of  mixtures  of  systematic  and  random  error  aggre¬ 
gates.  These  would  degrade  the  estimation  of  post-mission  corrections.  The' 
may  be  obtained  in  the  following  way:  let  E0  +  and  %e  +  3Se  be  observab! 
initial  and  terminal  vertical  deflection  components,  D.  ■+■  A?e  the  GGI-deter-- 
mined  difference  between  these  deflections,  and  H0,  5E,j,  and  AEe  uncorrelat'  : 
random  errors.  The  problem  is  to  estimate  the  variable  correction  in  the  for:' 
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and  two  corresponding  additional  equations  with  5n(t)  and  5Ag(t)  to  be  esti¬ 
mated.  The  first  of  three  covariance  equations  reads  then,  with  bars  indicat 
ing  covariances. 


— 2i  :  — 2  — : 

(  )  *  a  .  ,  2(  i")  *  A-l  in 

e  12 ,  e 


(23 


where  2  var  ig  has  been  omitted  because  of  the  relatively  high  vertical  grav¬ 
ity  component  measurement  accuracy.  As  an  alternative  post-mission  adjustment 
metnod,  empirical  corrections  applying  to  fixed  traverse  intervals  may  be 
calculated  under  determination  of  a  few  polynomial  constants  under  considera¬ 
tion  of  independently  provided  gravity  vector  information.  Metzger  and 
Jincitano  (1982)  have  published  preliminary,  non-optional  deflections  rms 
error  estimates  summarized  below  under  neglect  of  initial  and  terminal  astro- 
geodetic  deflection  errors. 


TABLE 

Sell-Estimated  Pore  Gradiometric  Deflection  RMS  Errors 
for  a  Tie  Point  Separation  of  80  Km  and  Vehicle  Speed  of  30to  150  Kmhr-' 


True*  Length  (Km) 

o 

20 
~»t  i 

ou 

8'.' 


llong  Track  Error 

0.000 
0.053 
0.06  3 
0.O53 
0.00  0 


(arcsec)  Cross 


Track  Error  (arcsec; 

0.000 

0.O91 

0.111 

0.091 

0.000 
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A  reduction  of  GGI-generated  systematic  and  random  errors  can  be  achieved  by 
intentional  gradiometer  platform  motion,  primarily  about  the  vertical  spin 
axis.  This  results  in  modulated  signal  outputs.  White  (1980)  established 
that  the  rotation  or  carousel  frequency  wc  >  ir/^V,  where  V  designates  vehicle 
velocity.  The  minimum  carousel  frequency  for  a  land  vehicle  with  speed  10 
msec  would  thus  be  3.6  cycles  hr”  .  For  plane  velocities,  the  minimum 
requirement  becomes  prohibitive.  For  aircraft  applications  under  coverage  of 
areas  of  the  order  (300  km)  ,  parallel  and  cross  traverses  separated  by  about 
5  km  and  multiple  tie  points  are  a  prerequisite  for  high  accuracy  gravity 
vector  determination.  To  minimize  the  effect  of  hight  attenuation,  the 
aircraft  altitude  should  not  exceed  a  height  of  500  m  above  terrain  level. 

The  downward  analytical  continuation  error  of  gravity  vector  components 
amounts  to  3Z  and  is  expected  to  be  somewhat  greater  over  strongly  mountainous 
terrain. 

3.  COMBINED  USE  OF  GRADIOMETER  AND  INERTIAL  SURVEY  DATA 

The  availability  of  inertial  data  at  vehicle  stops  during  terrestrial  surveys 
permits  simultaneous  astrogeodetic-gradiometric  and  astrogeodetic-inertial 
gravity  vector  determination.  Despite  an  inherent  greater  accuracy  of  the 
gradioraetric  method,  the  inertial  method  can  profitably  serve  control  pur¬ 
poses.  In  the  case  of  acceptable  agreement,  weighted  averages  may  be  computed 
under  use  of  covariance  analysis.  Because  of  Che  inclusion  of  linear  combina¬ 
tions  and  correlations  of  astrogeodetic  deflections,  a  weighted  mean  of  the 
rorm  var^  (varj  *  va^)”1  x,  +  var^  (varj  +  va^)”^  cannot  be  employed, 
i.e.,  che  second,  ''inertial*  weight  factor  becomes  relatively  smaller  through 
covariance  analysis.  Further,  the  weights  differ  from  point  to  point.  On  Ch. 
other  hand,  astrogeodetic  deflections  with  rms  errors  of  0.25  arcsec  make  the 
inertial  multiple  track  method  attractive  since  the  variable  astrogeodetic- 
gradioraecer  and  astrogeodetic-inertial  error  variances  would  not  differ  appre¬ 
ciably. 

6.  APPLICATIONS 

The  astrogeodetic-inertial  and  astrogeodetic-gradiometric  methods  allow  the 
determination  of  the  detailed  structure  of  the  gravity  field  in  continental 
and  adjacent  areas  with  the  former  method  restricted  to  semi-flat  terrestrial 
surveys.  Satellite-based  gravity  gradiometers  with  spatial  deflection  control 
data  over  land  areas  could  provide  the  means  for  the  computation  of  gravity 
potential  coefficients  of  about  degree  and  order  90  in  rather  inaccessible 
areas.  Gravity  tensor  information  in  selected  regions  would  facilitate  fine 
structure  analyses.  Large  scale  or  small  scale  utilization  of  the  available 
technology  would  result  in  improved  space  vehicle  trajectories  and  orbits, 
global  positioning  and  flight  navigation,  gravity  programmed  inertial  posi¬ 
tioning  systems,  subterraneous  mass  detection,  and  scientific  investigations 
inc Lading  cests  of  Newton's  mass  attraction  law  by  means  of  the  superconduct¬ 
ing  gravity  gradiometer. 

7.  CONCLUSIONS 

oitton  and  Honeywell  test  data  produced  by  essentially  first  generation  iner¬ 
tial  sytems,  existing  gyros,  accelerometers  and  velocity  quantizers  having 
onsiierably  higner  performance  characteristics  than  presently  installed 
.  istruments,  and  advanced  optimization  methods  indicate  potential  for  gravity 
•vitir  component  determination  with  rms  accuracies  of  about  3.3  arcsec  and  0.3 
mgil  <r  somewhat  better  in  the  context  of  area  adjustments  under  use  of 
:u.ti,ue  istr  jgeodet  u  deflections  with  sta  t  i  s  t  ioa  1  ly  independent  errors.  The 
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astrogeodetic-inertial  method  Is  In  general  sufficiently  accurate  in  semi-flat 
terrain.  The  Wiener-Kolmogorov  horizontal  channel  deflection  determination 
method  is  optimal  for  single  and  multiple  track  data  and  may  be  slightly 
modified  or  perfected  for  most  profitable  application.  The  Bell  gravity 
gradloaeter  is  efficient  and  indispensable  for  large  area  aircraft  surveys  and 
over  strongly  mountainous  terrain.  Optimal  post-mission  adjustment  methods 
have  to  be  perfected  in  connection  therewith.  Combined  astrogeodetic- 
gradiometric  and  astrogeodetic-inertial  terrestrial  surveying  is  feasible  and 
advantageous  for  control  purposes.  Appropriate  weighting  factor  computations 
require  covariance  analyses.  The  identified  advanced  systems  and  mathematical 
methods  increase  significantly  the  state-of-the-art  in  physical  geodesy  and 
surveying  and  permit  multiple  applications. 
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